% 2-D IAEA reactor benchmark, est. k = 1.0296

it = 1;
buck = 0.8e-4;
dat = [ 1.500   0.030   0.000  1.000  0.000  0.000  
        0.400   0.080   0.135  0.000  0.020  0.000
        1.500   0.030   0.000  1.000  0.000  0.000
        0.400   0.085   0.135  0.000  0.020  0.000
        1.500   0.030   0.000  1.000  0.000  0.000  
        0.400   0.130   0.135  0.000  0.020  0.000 
        2.000   0.040   0.000  0.000  0.000  0.000
        0.300   0.010   0.000  0.000  0.040  0.000
      999.001   0.000   0.000  0.000  0.000  0.000
      999.001   0.000   0.000  0.000  0.000  0.000 ];     
dat(:,2)=dat(:,2)+buck*dat(:,1);    

numg = 2;  numm = 5;
xcm  = [  0   10   30   50   70   90  110  130  150  170];
xfm  = [     10   20   20   20   20  20   20   20   20 ]/1;
ycm  = [   0   10   30   50   70   90  110  130  150  170  ];
yfm  = [     10   20   20   20   20  20   20   20   20 ]/1;

% mt   = [  3  2  2  2  3  2  2  1  4  
%           2  2  2  2  2  2  2  1  4
%           2  2  2  2  2  2  1  1  4
%           2  2  2  2  2  2  1  4  4 
%           3  2  2  2  3  1  1  4  5
%           2  2  2  2  1  1  4  4  5
%           2  2  1  1  1  4  4  5  5
%           1  1  1  4  4  4  5  5  5
%           4  4  4  4  5  5  5  5  5 ];
      
mt   = [  3  2  2  2  3  2  2  1  4  
          2  2  2  2  2  2  2  1  4
          2  2  2  2  2  2  1  1  4
          2  2  2  2  2  2  1  4  4 
          3  2  2  2  3  1  1  4  4
          2  2  2  2  1  1  4  4  4
          2  2  1  1  1  4  4  4  4
          1  1  1  4  4  4  4  4  4
          4  4  4  4  4  4  4  4  4 ];      

BCL = 1;
BCR = 0;
BCB = 1;
BCT = 0;        

src(:,:,1) = [ 0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0 
               0  0  0  0  0 ];            
src(:,:,2) = [ 0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0 
               0  0  0  0  0 ];           
src(:,:,3) = [ 0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0 
               0  0  0  0  0 ]; 

IBSL = [0 0 0 0 0
        0 0 0 0 0]';
IBSR = [0 0 0 0 0
        0 0 0 0 0]';
IBSB = [0 0 0 0 0
        0 0 0 0 0]';
IBST = [0 0 0 0 0
        0 0 0 0 0]';